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Abstract 

We present a general hydrodynamic theory of transport in the vicinity of superffuid-insulator 
transitions in two spatial dimensions described by "Lorentz" -invariant quantum critical points. We 
allow for a weak impurity scattering rate, a magnetic field B, and a deviation in the density, p, from 
that of the insulator. We show that the frequency-dependent thermal and electric linear response 
functions, including the Nernst coefficient, are fully determined by a single transport coefficient 
(a universal electrical conductivity), the impurity scattering rate, and a few thermodynamic state 
variables. With reasonable estimates for the parameters, our results predict a magnetic field 
and temperature dependence of the Nernst signal which resembles measurements in the cuprates, 
including the overall magnitude. Our theory predicts a "hydrodynamic cyclotron mode" which 
could be observable in ultrapure samples. We also present exact results for the zero frequency 
transport co-efficients of a super symmetric conformal field theory (CFT), which is solvable by 
the AdS/CFT correspondence. This correspondence maps the p and B perturbations of the 2+1 
dimensional CFT to electric and magnetic charges of a black hole in the 3+1 dimensional anti-de 
Sitter space. These exact results are found to be in full agreement with the general predictions of 
our hydrodynamic analysis in the appropriate limiting regime. The mapping of the hydrodynamic 
and AdS/CFT results under particle- vortex duality is also described. 
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I. INTRODUCTION 



A key indication that the normal state of the cuprate superconductors is aberrant came 
from the pioneering measurements of the Nernst effect by Ong and collaborators 
Also of interest here are measurements of the Nernst effect in Nbo.15Sio.s5 films by Behnia 
and collaborators 0, 0]. The Nernst co-efficient measures the transverse voltage arising in 
response to an applied thermal gradient in the presence of a magnetic field. The response 
of Fermi liquids is associated with a weak particle- hole asymmetry in the spectrum of the 
fermionic excitations near the Fermi level [7|. The large observed response, and its striking 
and unexpected dependence on the magnetic field, temperature, and carrier concentration 
indicated that an explanation starting from a metallic Fermi liquid state could not be tenable. 
Instead, Ong and collaborators argued that their observations called for a description in 
terms of a liquid of quantized vortices and anti-vortices in the superconducting order (and 
its precursors) at low temperatures. 

A complete theory of the dynamics of the vortex liquid state is so far lacking. Ussishkin 
et al. 8] used a classical Gaussian theory of superconducting fluctuations, and Mukerjee and 
Huse 9] extended this to a time-dependent Ginzburg-Landau model. Podolsky et al. [h]] 
applied classical Langevin equations to a model of phase variables residing on the sites of a 



hypothetical lattice. Anderson has taken a speculative view of the vortex liquid, arguing 
against the conventional Debye-screening of vortex interactions. While some experimental 
trends are successfully described by Refs. , [h], it would be useful to have a kinematic 
approach which is also able to include quantum effects, and extends across the superfluid- 
insulator transition. Quantum effects which surely play an more important role at lower 
temperatures, especially in the underdoped region. Indeed, it is the equal importance of 
thermal and quantum fluctuations which underlies the difficulty in describing this vortex 
liquid. 

As in the recent work by Bhaseen, Green, and Sondhi 12|, this paper will advocate an 
approach departing from the quantum critical region of a zero temperature (T) quantum 
phase transition between a superconductor and an insulator. This is the region where 
the primary perturbation from the physics of the T = quantum critical point is the 
temperature. The single energy scale, ksT, then determines observable properties, including 
the values of diffusion constants and transport co-efficients, in a manner that has been 
discussed at length elsewhere [l^, Q. The electrical conductivity of this quantum critical 
system, which we denote ctq, will play a prominent role in our results. In 2+1 dimensions, 
near quantum critical points which obey hyperscaling properties, this conductivity is given 
by EQEl 

4e 2 

a Q = (1-1) 

where in the quantum critical region $ CT is a universal dimensionless number dependent only 
upon the universality class of the critical point. 
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FIG. 1: Zero temperature (T = 0), zero field (B = 0) phase diagram in the vicinity of the 
quantum critical point described by the CFT, represented by the filled circle. The coupling g 
represents a parameter which tunes between a superfluid and a Mott insulator which is at a 
density commensurate with the underlying lattice. The chemical potential p introduces variations 
in the density and p is difference in the density of pairs of holes in the superfluid from that in the 
Mott insulator. The thin dashed lines are contours of constant p. In the application to the cuprate 
superconductors, the Mott insulator with p = could be, e.g., an insulating state at hole density 
5i = 1/8 in a generalized phase diagram; then p = [8 — 5j)/(2a 2 ), where a is the lattice spacing. 
The thick dotted line represents a possible trajectory of a particular compound as its hole density 
is decreased; note that the ground state is always a superconductor along this trajectory, even at 
5 = 1/8 (although there will be a dip in T c near 5 = 1/8 as is also clear from Fig. [2]). Note that the 
parent Mott insulator with zero hole density is not shown above. This paper will describe electrical 
and thermal transport in the above phase diagram perturbed by an applied magnetic field B and 
a small density of impurities. 



The discussion so far applies, strictly speaking, only to systems which are exactly at the 
commensurate density for which a gapped Mott insulator can form. The cuprates, and other 
experimental systems, are not generically at these special densities, and so it is crucial to 
develop a theory that is applicable at generic densities. Such a theory will emerge as a 
special case of our more general results below. We allow the density to take values p by 
applying a chemical potential shown in Figs. [1] and EJ 

We emphasize that p measures the deviation in particle number density from the density 
of the commensurate insulator 19|, |20|; so p can be positive or negative, and we will see 
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FIG. 2: Nonzero temperature (T) phase diagram at B = along three vertical cuts {i.e. fixed g) 
in Fig. [TJ The lines indicate Kosterlitz-Thouless phase transitions at T = Tkt associated with the 
loss of super fluid order as a function of p for different values of g. At g = g c , Tkt/\p\ is a universal 
number determined by the CFT at g = g c , p = [3]. This paper describes transport properties 
in the non-super fluid region above Tkt, in the presence of an applied magnetic field B and a 
small impurity scattering Ti mp . The results of the supersymmetric CFT solvable by AdS/CFT in 
Section [V] are limited to g = g c , but allow arbitrary variations in p and B away from quantum 
criticality as long there is no phase transition into a superfluid (or other) state. 

that the sign of p determines the sign of the Hall resistivity and other transport coefficients. 
Also, purely as a choice of convention, we will measure p in terms of density of pairs of 
holes; this choice does not imply the degrees of freedom of the underlying theory are Cooper 
pairs, although this is the case in the simplest model. We will use general hydrodynamic 
arguments (specialized to "relativistic" quantum critical points) to show that the frequency 
(u) dependent conductivity, o xx at a generic density p is given by 



where the system is characterized by the thermodynamic state variables e, the energy density, 
and P, the pressure-we will specify their temperature dependencies shortly in Eq. (11.81) . The 
factor of (2e) 2 is a consequence of our choice for the normalization of p\ note that product 
2ep measures the net charge density, and so is independent of this convention. We assume 
there is a dilute concentration of impurities which relax the conserved momentum 1 , and lead 
to the scattering rate l/r imp — the temperature dependence of r imp is specified in Eq. (11.111) . 



Umklapp scattering can also lead to momentum relaxation, but this is exponentially small at low T. 



a 



XX 
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The parameter v is a velocity characteristic of the quantum critical point which is assumed 
to have dynamic critical exponent z — 1. Finally, the crucial parameter ctq, is the same 
quantum conductance which appeared in Eq. fll.ip . However, as one moves away from the 
critical coupling g = g c and adds a non-zero p in Fig. [21 $ CT will acquire a dependence upon 
the ratios (g — g c ) /T 1 ^ and p/T which can be included unchanged in our results below (here 
v is the usual correlation length exponent). 

It is perhaps helpful to note here the "non-relativistic" limit of Eq. (II. 2p . which does not 
constitute the regime of primary interest of this paper. In this limit e + P ~ \p\mv 2 (where 
m is the mass of the particles), and then the second term takes the form of the conventional 
Drude result. 

To develop a theory for the Nernst effect, we need to apply a magnetic field B to the 
system described so far. A central result of this paper is that for not too strong B fields, 
the Nernst response, and a set of related thermoelectric transport co-efficients, are com- 
pletely determined by the thermodynamic variables and impurity scattering rate appearing 
in Eq. (11.21) and the single universal transport co-efficient oq. In particular, no additional 
transport co-efficients are needed. Thus there are a large number of Wiedemann- Franz-like 
relations which relate all the thermoelectric response functions to the regular part of the 
electrical conductivity in zero field, o~q. We will also determine the frequency dependence 



left expressions are given below. 

121 ] only considered a non-zero B, with p = l/r imp 



of these transport co-efficients; exp 

In their work, Bhaseen et al. 
0. Their primary new result concerned the longitudinal thermal conductivity, k xx at zero 
frequency. Our result for k xx (uj = 0) is consistent with theirs, and further, we show further 
that it is related to <jq by a Wiedemann- Franz like identity (Eq. (I1.28P below). However, 
remarkably, unlike the conventional identity which specifies the ratio of k xx to the electrical 
conductivity, our identity specifies the ratio of k xx to the electrical resistivity. This suggests 
a physical picture of transport currents carried by vortices rather than particles, at least 
when the perturbation associated with B is larger than that associated with p. 



A. Characterization of systems under consideration 

Let us now specify the class of theories to which our results apply. Current theories of 
the superfluid-insulator transition in non-random systems in 2+1 dimensions are described 
by quantum field theories which are Lorentz invariant, and are therefore conformal field 
theories (CFTs). Consequently, we will mainly restrict our attention here to T > "quantum 
critical" phases of CFTs, and the general structure of their response to a non-zero p (which 
is not restricted to be small) and small B and l/r imp . We expect that many of our results, 
and especially the magnetohydrodynamic analysis in Section IHIt can be generalized to a 
wider class of systems, but we will not discuss such a generalization here. 

Specific examples of CFTs to which our results apply are: 
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(i) The superfluid-insulator transition of the boson Hubbard model on a two-dimensional 
lattice with a density of an integer number of bosons per unit cell. The bosons carry charge 
±2e because they are caricatures of Cooper pairs. The critical point is described [21| by the 
Wilson-Fisher fixed point of the \ifj\ A field theory of a complex scalar ip (representing the 
boson annihilation operator), see Eq. (II. 3p below. This field theory also has a dual repre- 



sentation [22|, 123( in terms of a vortex field <p coupled to an emergent U(l) gauge field. Our 
results apply equally to both representations, and the observable properties do not depend, 
naturally, on whether the particle or vortex representation is used to describe the CFT. 
(ii) The superfluid-insulator transition of the boson Hubbard model on a two-dimensional 
lattice with a mean density of a rational number, p/q (with p, q, co-prime integers), of 
bosons per unit cell. A 'deconfined' critical point [24j] is then possible JjJ, |25j between the 
superfluid and an insulator with valence-bond-solid order and is described by the theory of 
q flavors of vortex fields (ft coupled to an emergent U(l) gauge field. This field theory can 
also be 'undualized' to a "quiver gauge theory" of fractionalized bosons with charge ±2e/q 



19|. 



[iii) Electronic models with a <i-wave superconducting ground state can also undergo de 



confined phase transitions to insulating states with valence-bond-solid order |2_6|, 123] • The 
CFTs of these transitions have Dirac fermion degrees of freedom, representing the gapless, 
Bogoliubov quasiparticle excitations of the <i-wave superconductor, in addition to the mul- 
tiple vortex and gauge fields found above in (ii). 

(iv) Yang-Mills gauge theories with a SU(iV) gauge group and Af = 8 supersymmetry. These 
are attracted in the infrared to a superconformal field theory (SCFT) which is solvable in 



the 

I 



arge N limit via the AdS / CFT mapping. This solution has been used in previous work 



281 ] to obtain the collisionless-to-hydrodynamic crossover in the transport of a conserved 



SO (8) R-charge, as well as an exact value for <& a . Here we will examine, as in other recent 
work the deformation of the SCFT by a non-zero B and p. The B field and density 
p are both associated with a U(l) subgroup of the SO(8) R charge. After the AdS/CFT 
mapping, B and p correspond to the magnetic and the electric charge of a black hole in 
AdS space. We will present exact results for the conserved current correlators of the dyonic 
black hole in Section |Vj which allows us to obtain corresponding exact results for the Nernst 
and related thermoelectric responses of the SCFT. In the appropriate hydrodynamic limit, 
these results are found to be in full agreement with the more general magnetohydrodynamic 
analysis in Section IHIl Additional comparisons between the hydrodynamic and AdS / CFT 
results appear in a separate paper [30 ]. 

It is worth reiterating that not all of the above CFTs are purely bosonic, and the examples 
in (iii) and (it;) contain fermionic degrees of freedom. Furthermore, in cases (ii) and (iii), 
the bosonic degrees of freedom of the CFT are not Cooper pairs, but fractions of a Cooper 
pair with charges determined by the density of the Mott insulator. 
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B. A simple model 



Before presenting our main results, it is useful to establish notation by explicitly writing 
down the simplest of the CFTs listed above. This is the l^l 4 field theory for bosons with 
charge ±2e and action 



S 



d 2 rdr 



V - i^-A ] 

he j 



-<#i 2 +§i# 



(1.3) 



where r = (x, y) is a 2-dimensional spatial co-ordinate, r is imaginary time, g is the coupling 
which tunes the system from the superfluid to the insulator (see Figs. ITO) . and the quartic 
coupling u is attracted to the Wilson-Fisher fixed point value in the infrared limit associated 
with the CFT. The electromagnetic gauge potential is non- fluctuating (and is not to be 
confused with the emergent U(l) gauge field of the vortex CFTs noted above). Its time- 
component takes an imaginary value (in imaginary time) which determines the chemical 
potential 

i2eA T = p, (1.4) 



while the spatial components take r-independent values so that 



V x A = B, 



;i.5) 



with a spacetime-independent magnetic field B. The density, p, is defined, as usual, by the 
derivative of the partition function with respect to the chemical potential 



P 



k B T 
hV 



ds_ 

dp 



;i.6) 



where V is the volume of the system. We reiterate that p measures the difference in the 
density from that of the commensurate, T = 0, insulating state, and not the total density. 
Also, p is a charge density in the sense that it measures the number density of particles 
minus the number density of anti-particles. 

Another parameter above which will be important for experimental comparisons is the 
velocity v. Note that it plays the role of the velocity of "light" in the "relativistic" CFT. It is 
determined here by the parameters of the underlying boson Hubbard (or other microscopic) 
model whose superfluid-insulator transition is described by the above CFT. It is important 
to distinguish v from the velocity c, which is the actual velocity of light. Here c merely 
plays the role of a coupling constant which relates the value of B to physical CGS units, and 
is not a velocity associated with the dynamics of the physical model under consideration. 
Because « < c, we can neglect the actual relativistic corrections associated with the physical 
quantum fluctuations of the photon field A^. 

With the definition of v at hand, we can now begin comparing the various energy scales 
which characterize the system. The largest energy scales which characterize the deviation 
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from the T = quantum critical point are k B T, an energy scale mo ~ \g — gc\ u associated 
with the deviation from critical coupling, and the chemical potential p. We will generally 
assume that k B T is the largest of these scales; our results allow m and p to be of order k B T, 
but not too much larger. For the solvable SCFT theories considered in Section |V] the energy 
scales associated with p and B will not be restricted to small values. However, for the more 
general analysis in the remainder of the paper which applies also to non-supersymmetric 
CFTs, we will assume that the perturbation due to B is small, and in particular, 

hv^2eB/(hc) < k B T. (1.7) 

Some thermodynamic state variables will also appear in our transport result. Their tem- 
perature dependencies obey scaling forms similar to those computed earlier for the present 
theory at p = in Ref. Ell, and for p ^ in Ref. 32. In particular, we will need results for 



the energy density, e, and the pressure, P which obey 



2 



/ 1^1 <l> 

hv 



2 



P = k B T J *p, (1-8) 

where, as in Eq. (11.11) . $ £j p are dimensionless universal numbers which depend on the ratios 



(g — g c ) I \k B T) 1 / v and p/(k B T) 3l|, [32J. The dependence on B is not difficult to account 
for, but will be subdominant, and non-singular, under the condition in Eq. (11. 7ft . 

The final parameter to introduce in our theory of the Nernst effect and other thermoelec- 
tric response functions is the momentum relaxation rate l/r imp . The theory defined so far 
conserves total momentum, and thus, such relaxation requires an additional perturbation. 
We assume that the relaxation is caused by a weak random impurity potential V(r) which 
couples to the most relevant perturbation allowed by symmetry at the CFT fixed point. For 
the present theory, this is the operator \ip\ 2 , and therefore the impurity action is 

S imp = j dr j d 2 rV(r)\ij(r,T)\ 2 . (1.9) 

We will take a "quenched" average over the ensemble of impurity potentials which obey 



V(r) = ; V(r)V(r>)=V 2 ip S 2 (r-r'), (1.10) 



and work to order V| mp . Note that total energy and charge are conserved in the presence 



of V(r), and momentum is the only conserved quantity which will relax. We estimate the 
momentum relaxation rate in Section [IV] and find 

— -vl^- 2 /". (1.11) 



7~imp 



The condition for this scattering to be small is 

h/r imp < k B T. (1.12) 
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The present model has v ~ 2/3, and so l/r imp depends upon temperature only very weakly. 
Indeed, all the CFTs noted earlier are expected to have a similar value for v. It is therefore 
a reasonable first approximation to treat l/r imp as a temperature-independent constant. We 
will also ignore the dependence of l/r imp on B and p, under the condition in Eq. (11.71) . 



C. Summary of results for the thermo-electric response 

We finally turn to a statement of our main results for the transport co-efficients. We are 
interested in the response of the electrical current J and the heat current Q to an applied 
electric field E and a temperature gradient VT. The precise definitions of J and Q appear 
in the contexts of the models studied in the body of the paper. The electric field can be 
applied by allowing for a weak spatial dependence in the chemical potential \i (which is then, 
formally, the electrochemical potential) with 2eE = — V/i, while the temperature gradient 
describes a similar weak spatial dependence in T. The transport co-efficients are defined by 
the relation 

' J \ _ I a a \ I E 
Q J \Ta X [ -VT 



(1.13) 



where a, a and k are 2x2 matrices acting on the spatial indices x, y. Rotational invariance 
in the plane imposes the form 

a = a xx 1 + a xy e, (1.14) 

where 1 is the identity, and e is the antisymmetric tensor e xy = —e yx = 1. o xx and o xy 
describe the longitudinal and Hall conductivity, respectively. An analogous form holds 
for the thermoelectric conductivity a (which determines the Peltier, Seebeck, and Nernst 
effects), as well as for the matrix 7t which governs thermal transport in the absence of electric 
fields. The latter applies to samples connected to conducting leads, allowing for a stationary 
current flow. In contrast, the thermal conductivity, k, is defined as the heat current response 
to — VT in the absence of an electric current (electrically isolated boundaries). It is given 
by 

k = K — Tacr^a. (1-15) 

Finally, the Nernst response is defined as the electric field induced by a thermal gradient in 
the absence of an electric current, and is given in linear response by the relation E = — i?VT, 
with 

$ = -a- 1 a. (1.16) 

The Nernst signal is the transverse response, = fiyx- The Nernst co-efficient is usually 
defined as v = e^/B, which tends to become field independent at small B. The Nernst 
signal is expected to be positive if it is due to driven vortices, while it is generally negative 



if it arises from quasiparticle excitations 33| 
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We now present our main results for the transport co-efficients. For the computations 
using AdS/CFT applied to the SYM theory in Section [V] results can be obtained for general 
external frequency, uj. However, our more general hydrodynamic results apply only for 
hjj <C ksT, and this condition is assumed in the remainder of this section. We begin by 
presenting our complete result for the frequency dependence of the longitudinal electrical 
conductivity (whose B — > limit was already quoted in Eq. (II. 2p ): 



&xx — °"<2 



(u + if r imp ) (u + ij + ia;g/7 + i/r imp ) 
(lj + 27 + i/r imp ) 2 - u 2 



(1.17) 



The overall scale is set by the quantum conductance oq introduced in Eq. ( 11. 2p . and the 
remainder depends upon two important frequency scales which will appear throughout our 
analysis. These frequencies are 

2eBpv 2 



a ' s *fp>- (L18) 



and 



V 



^WTpy (L19) 

We identify the first frequency, u c , as a cyclotron frequency. This seems a natural inter- 
pretation in view of the damped resonance present in the denominator of Eq. (I1.17p . Note 
that in the non- relativistic limit where e + P ~ \p\mv 2 , u c reduces to the familiar result 
u c = 2eB/(mc). For relativistic particles the cyclotron frequency decreases with the energy 
E as uj c (E) = 2eBc/E. In the present context where v plays the role of the velocity of 
light, this is modified to u c (E) = 2eB/c-v 2 /E. The hydrodynamic expression f 1 1 . 1 8 j) can be 
regarded as a thermal average over cyclotron frequencies uj c (E ~ T), while the proportion- 
ality to the charge density, u c ~ p, reflects the fact that particles and antiparticles circle in 
opposite senses. 

We can consider the cyclotron mode as arising either from the motion of particles and 
anti-particles, or from the motion of vortices and anti-vortices. In the latter interpretation, 
the roles of B and p are interchanged, whilst the expression for u c remains invariant. We 
will have more to say about this 'dual' interpretation in the body of the paper, and further 
results appear in a separate paper j3o| . 

The second frequency, 7, is the damping frequency of the cyclotron mode of particles and 
anti-particles. Note that this damping is present even in the absence of external impurities, 
and is a consequence of collisions between particles and anti-particles which are executing 
cyclotron orbits in opposite directions. This should be contrasted from the behavior of a 
Galilean-invariant system (i.e., a system with no anti-particle excitations) for which Kohn's 



theorem 34J guarantees an infinitely sharp cyclotron mode. The sharpness of the cyclotron 



resonance is determined by the ratio 



7 B 

- = — , (1-20) 

u c (p p 
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which up to the factor & a equals the number of flux quanta, 

he 

<P„ = Ye , (1.21) 

applied per charge 2e. 

We will see later that a different frequency plays the role of the damping of the cyclotron 
mode when it is interpreted as due to the motion of vortices and anti-vortices. In that case 
the damping frequency is 

u 2 c 4e 2 pV 

There is an obvious 'dual' structure apparent upon comparing Eq. (11.191) and (11.221) . which 
we will discuss in more detail. Note that the cyclotron resonance will be visible only in 
ultrapure samples where l/ri mp <C lo c . In this case, the cyclotron resonance is sharp in the 
thermoelectric response associated with particle transport if 7/cj c = B/(<pop) <C 1, while the 
same condition implies a washed out resonance in the dual response functions associated 
with vortices. In the opposite regime, 7/o; c ^> 1, the vortex response, and in particular the 
Nernst effect, should exhibit a sharp cyclotron resonance in ultrapure samples. 

Another notable feature of Eq. (11.171) is the singular nature of the limit associated with 
the small perturbations of the quantum critical region of Fig. [2J In particular, note that for 
the d.c. conductivity at to = 0, the value of a xx depends upon the order of limits of p — ► 0, 
B — > 0, and l/ri mp — > 0. This singular limit reflects the fact that the low frequency transport 
studied previously at p = B = in Refs. has ballistic energy propagation and an 

infinite thermal conductivity. For non-zero p or B, the energy and number currents can mix 
with each other, leading to a finite thermal conductivity and an order unity correction to 
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a xx , as anticipated for the case p = in Ref. 

A useful property of Eq. (I1.17P is that it depends only upon the combination to + i/r imp . 
This is actually a property of the long distance limit of the hydrodynamic equations presented 
in Section UTT| and is obeyed by all the transport co-efficients. The remainder of this section 
will therefore present results only in the d.c. limit, while the u dependence can be easily 
reconstructed by replacing l/r imp — > l/r imp — iuj (as long as tvjj <C ksT). 



D. Nernst effect 

Our central result for the Nernst signal is 

6 AT = $yx = 



k B \ fe + P 



imp 



2e J \ k B Tp 
1 fk B \ ( e + P 



{ulh + l/r imp ) 2 + ^1 



imp 



1.23) 

;i.24) 



$ a \2e) \k B TB/<p ) L(^/ 7 + l/r imp ) 2 + W 2 

where Q a is the dimensionless universal number appearing in the expression for the con- 
ductivity <jq in Eq. (jl.ip . We have expressed the Nernst signal in terms of its quantized 
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unit, 

= 43.086 fiV/K, (1.25) 

times dimensionless ratios in the various brackets. We can use the relation e+P ~ Ts, where 
s is the entropy density (see Eq. (13.151) ). valid for small p, p, to identify the factor in the 
second brackets as approximately the entropy per particle in Eq. ( 11.231) . and as the entropy 
per vortex in Eq. ffl~2l) . The combination of Eqs. ([HI]), (Oil . (jTTTj) and 01. 2301. 241) now 
implies an interesting and non-trivial dependence of the Nernst signal on B and T. Those 
should be observable in experimental regimes where the entire thermoelectric response is 
dominated by critical superconducting fluctuations, as will be discussed further in Section HT1 



E. Other thermo-electric transport coefficients 



We conclude this introductory section by mentioning two other results for transport 
co-efficients whose limiting forms can be compared with earlier computations. For the 
transverse thermo-electric conductivity we obtain 



a 



XL) 



2ek B \ ( s/k B 
h J \B/<p 



1 2 + uj 2 c + 1 /r imp (l-pp/(Ts))' 



(7 + l/r imp ) 2 + u, 



1.26) 



While in most recent experiments, the electric conductivity a receives the largest contribu- 
tion from non-critical carriers, the thermoelectric conductivity is dominated by supercon- 
ducting fluctuations, even far above T c . It is thus the main quantity to be compared with 
recent experimental observations in Section [III As with earlier results, a xy has been written 
in terms of the quantum unit of the thermoelectric co-efficient p, [lo|, 2eks/h = 6.7 nA/K, 
times dimensionless ratios. In the absence of impurity scattering, l/r imp — > 0, the factor 
in the square brackets is unity, and we have ot xy = sc/B, a result obtained long ago for 
non-interacting fermions 
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MM and later derived by Cooper et al |43) for interacting 
fermions, and by Bhaseen et al. [12J for the superfluid-insulator transition. 



For the longitudinal thermal conductivity we obtain 



M 






1 









p 



2 r 



(a; c 2 / 7 )K 2 /7+l/ri 



imp I 



e + P 
k B TB/& 



(a; c 77+l/r imp ) 2 + a;, 



(w 2 /7 + l/r imp ) 2 + ^ CJ 



(1.27) 
(1.28) 

where now k B T/h is the quantum unit of thermal conductance. In the limit l/r imp — > 
and B — > 0, the factor within the square brackets in Eq. (I1.27P reduces to unity. The 
resulting expression for k xx relates it to <jq in a Wiedemann-Franz-like relation, as has been 
noted by Landau and Lifshitz j^fj (and elaborated on recently in Ref. 0). This relation 
suggests a physical picture of transport due to particles/anti-particles carrying charges ±2e 
and entropy per particle s/p. 
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In the complementary limit of l/T- imp — > and p — > the factor within the square brackets 
in Eq. (11.281) reduces to unity. Now k xx is proportional to the resistivity 1/ctq, indicating 
a picture of transport due to vortices of net density B/<f) . The value of _k^ x has the same 
dependence upon all parameters as that obtained by Bhaseen et al. [12]. We can also 
com par e the value of the numerical prefactor. For <E> CT we use the value 1.037/e 2 obtained in 
Ref. [l7| in the e-expansion (e = 3 — d where d is the spatial dimension), which is the same 
expansion by Bhaseen et al. [12j. It is also easy to compute the value of $ e + $p in the same 
expansion: to the leading order needed, these are just given by the values for free, massless, 
relativistic bosons in d = 3, which yields $ e + <3>p = 4jr 2 /45 + 0(3 — d). Using these values we 
obtain the same result for Eq. (24) of Ref. [12|, with their dimensionless parameter 

g = 4.66 (not to be confused with our coupling g). This is to be compared with their value 
g = 5.55. The origin of this numerical discrepancy is not clear to us. We believe Eq. dl . 28[) 
is an exact identity in d = 2, but it is possible that it is modified when d is close to 3. 

The outline of the paper is as follows. Section [TT] will compare the result for the Nernst 
effect and the thermoelectric response Eq. (11.26)) with experiments on the cuprate supercon- 
ductors and on Nbo.15Sio.85 films. Section [Till will present the derivation of these results using 
a hydrodynamic analysis of transport near a generic, 2+1 dimensional, "relativistic" quan- 
tum critical point perturbed by a chemical potential, a magnetic field, and weak impurity 
scattering. An estimation of the impurity scattering rate appears in Section [IVl Section fVl 
will describe the exact solution for transport near a supersymmetric quantum critical point, 
perturbed by a chemical potential and a magnetic field, which is solvable by the AdS/CFT 
mapping to the physics of a dyonic black hole in 3+1 spacetime dimensions. Some technical 
details appear in the appendices. 



II. COMPARISON WITH EXPERIMENTS 

Our main results for the Nernst signal have already been stated in Section II Dl In the 
following subsections, we will these results with recent observations in the cuprate super- 
conductors , and also briefly discuss experiments in Nbo.15Sio.s5 films Bfl. As 
mentioned before, in most of these experiments the electrical conductivity is dominated by 
non-critical fermionic contributions which are not captured by our relativistic hydrodynamic 
description. On the other hand, the transverse thermo-electric response a xy is expected to 
be predominantly due to superconducting fluctuations and the vortex liquid. In practice a xy 
is conveniently measure via the Nernst signal, using the relation a xy ~ o xx r d yx . The latter 
holds if the non-critical Hall conductivity is small, a xy <C cr xx , as is usually the case. 

It is convenient to perform the experimental comparisons by rescaling B and p so that 
they are both measured in units of (energy) 2 , 

B = B<f )Q /(hv) 2 ; p = p/(hv) 2 . (2.1) 

Further we observe that in typical experiments the flux per (excess) particle is very small, 
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B/p <C 1 and therefore 7/w c < 1. In this regime Eq. (11.26P simplifies to 



2ek B J ' B/(f)o P 1 + (n mp uj c y 



S / *; + p(W + B 2 7> 2 (2jrT lmI> /fi)2 



where in the second line we have assumed a fully relativistic regime with s ~ T 2 and 
e + P ~ T 3 , and pp <C sT. We recall that $ £ +p and $ CT are universal functions of p/T, and 
have an additional dependence on (g — g c )/T l l p ' . 



A. The cuprates 

Given the relative simplicity of our model of the cuprate superconductors, detailed quan- 
titative comparisons with the observations of Ref. Q are probably premature. In particular, 
we have omitted the influence of long-range Coulomb interactions, which modifies the spec- 
trum of boson density fluctuations, and likely leads to a superfluid-insulator quantum critical 



point which is not Lorentz invariant [35l. l36j] . Also, although Dirac fermion excitations are 
included in some of the CFTs mentioned above (corresponding to the nodal points of the 
c?-wave superconductor) other Fermi excitations associated with a Fermi surface may also be 
important, especially in the case of NbSi. Keeping these caveats in mind, it is nevertheless 
useful to examine the extent to which the present model can describe the observations. As 
we will now show, using physically reasonable values of the parameters in the theory, our re- 
sults describe the overall absolute magnitude of the observations, and numerous qualitative 
trends 0] remarkably well. 

We work here with a simple caricature of our predictions: We ignore the T and p 
dependence of the universal functions and simply treat them as constants, $ CT ~ 1.037, 
^e+p ~ ~ 3C(3)/7r ps 1.148. This is equivalent to assuming in Figs. [TJ2] that g = g c 
and p = p = for the purpose of evaluating these functions. It is not difficult to extend our 
theory to include the influence of these corrections to the leading quantum-critical behavior, 
but such a detailed analysis would not be commensurate with the other simplifications noted 
above. 

We notice that for small B, Eqs. (12. 2112. 3p predict a Nernst signal linear in B. At not 
too large temperatures, the second term in the numerator of (12. 3 j) can be neglected and the 
ratio a xy /B is seen to increase with decreasing temperature as 1/T 4 , 



ir {B -* 0) ~ Usv {—) W { ] 

Such a power law with exponent 4 was indeed observed over two orders of magnitude in signal 



strength in underdoped La 2 _sSr,5 Cu0 4 (LSCO, 5 < 0.12), cf., Fig. 4 in Ref. llOl. Assuming 
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20 40 60 80 100 
T[K] 



FIG. 3: Contour plot (with logarithmic spacing) of the thermoelectric conductivity a xy (Eq. 12. 3[) 
as a function of temperature T and magnetic field B, for parameters hv = 47 meV A, 5 — 5i = 0.025 
and Ti mp = 10 _12 s estimated for LSCO. In the ordered low temperature regime T < T c ~ 30K, 
Eq. (|2,3p will receive modifications. 



a typical doping 5 — 5j = —0.025 for underdoped LSCO with a lattice constant a = 3.78 
A, we obtain a constraint for TmnV 3 from fitting (12 .4p to the experimental value a xy /B = 
0.48/ (T/30K) 4 nA/KT per layer [lOj. Assuming a typical scattering time r imp ~ 10 _12 s, we 
obtain an estimate for the velocity hv ~ 47 meVA. These are reasonable parameter values, 



37. 



with the velocity v being of the order of the characteristic velocity found in Ref. 

The result of Eq. (12. 2p is plotted as function of both T and B in Fig. [31 This contour 
plot should be compared, e.g., with Fig. 13 in Ref. m in the underdoped regime, above the 
superconducting transition T c . 

Using the above parameter estimates we predict the cyclotron resonance 

u c = 6.2GHz- JL f^SV, (2.5) 



IT V T 

which, at T = 35K, is by a factor u c /ujc e ^ = 2m^v 2 ■ p/(s + P) ~ 0.035 smaller than the 
cyclotron frequency of free electrons. However, as mentioned before, this resonance can only 
be observed in ultrapure samples where l/r imp <^ u> c , which is clearly not the case in LSCO. 

Having estimated the velocity v and the scattering rate Ti mp , we can make a quantitative 
prediction for the Nernst signal in the vicinity of a quantum critical point where the entire 
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20 40 60 80 
T[K] 



FIG. 4: Contour plot (linear scale) of the Nernst signal e^v = $yx (Eq. 12. 6p close to a quantum 
critical point, as a function of temperature T and magnetic field B. The parameters are the same 
as for Fig. O The signal strength in the plot ranges up to 10^V/K. 

thermo-electric response is expected to be dominated by critical fluctuations. In this case, 
Eq. fl 1 . 2 3 [) can be cast into the form 



■& = — 

u yx I 



2e 



($ £ + ^ P y^B(k B Tf(h/(2TTT imp )) 



P 2 



$ CT ($ £ + $ P )(k B T)3(h/(2nT imp )) + $2p 2 5 



which is plotted in Fig. H] 



1. Hall resistance 



(2.6) 



Very recently, measurements of the Hall resistance in the high field normal state of 
YBa 2 Cu30 6 . 5 have been reported [38(. The focus of the authors was on magnetoresistance os- 
cillations; these oscillations are quantum interference effects which cannot be reproduced by 
the effective classical hydrodynamic models employed here (under the condition in Eq. (11.71) ). 
and so are beyond the scope of the present paper. However, the authors also reported a 
background Hall resistance which, surprisingly, was negative. The sample has hole density 
5 = 0.1. As argued in Section [Ij the density of mobile carriers, p, which appears in the 
hydrodynamic theory 19| (and which contributes a Magnus force on vortices 20j) is given 
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by the difference in density between the superconductor and the proximate Mott insulator. 
Using an insulator at 5j = 0.125, we have p = — 0.025/(2a 2 ). This negative value of p 
provides a very natural explanation of the observed negative Hall resistance. Also, we can 
predict that the Hall resistance should change sign as 5 is increased beyond 5p 

We can make a more quantitative comparison with experiments. In Eq. (I3.38P in Sec- 
tion IIII A[ we report the value of the Hall resistivity, p xy , and the zero frequency limit of 
that result is 



Pxy 



B 

2epc 



imp I 



(l/T- imp + ^/7) 2 + ^ 2 



(2.7) 



in the absence of impurity scattering (r; mp — > oo), this result was noted in Ref. [29j. Using 
the value of p noted above, at B = 60T, we determine that the prefactor of the square 
brackets in Eq. (12.71) is -2.2 kQ. For the factor within the square brackets, we assume the 
same parameters as found above for LSCO, and conclude that it is close to unity. This result 
is to be compared with the observed resistance per layer [38j] at this field of -3.9 kQ, which 
is quite reasonable agreement for this simple model. 



B. NbSi 



a™ oc — j OC ; ; — — -, (2.8) 



We also note the experiments on amorphous films of Nbo.15Sio.85 reported in Refs. 
The normal phase, T > T c , of these films exhibits a number of features that are consistent 
with our hydrodynamic results, when taken to their non-relativistic limit. In particular, 
ot xy as given in Eq. (12.21) displays a functional dependence on magnetic field akin to that 
reported in Ref. f]: 

B B 

OC 

C -4 + £ -4 l + {B/Bof 

with B = hc/el; 2 , which was interpreted as the physics being controlled by the shorter of 
the superconducting coherence length £ and the magnetic length £b — (hc/eB) 1 ^ 2 [6|. We 
mention that the low B data, i.e., the Nernst coefficient v = liniB^o a xy /Ba xx measured in 
Ref. H, was very successfully described by the theory of Gaussian fluctuations 0] . However, 
the crossover ( 12. 8ft and the high field behavior a xy ~ 1/B remained unexplained in earlier 
theories. Our magnetohydrodynamic approach may give a hint to the origin of the latter. 
We believe that the similarity of Eqs. (12. 2p and (12.81) is not a mere coincidence. Rather, it 
leads us to speculate that the scattering time r imp should be identified with 

rimp = = « (k F £) ~ r GL , (2.9) 



the Ginzburg-Landau life time of fluctuating Cooper pairs 39j. Here we have used the free 
electron value (non-relativistic limit) for the cyclotron frequency ujf^ = eB/m^c. Further, 
vp = hkp/m is the Fermi velocity, and we have used that kpi w 0(1) in the studied 
amorphous NbSi ^]. The estimate Ti mp ~ tgl suggests that the suppression of the Nernst 
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signal at high fields is due to the Cooper pairs starting to perform entire orbits over their 
life time. 

We may use the above guess of r imp to express the low field limit of a xy , as 

a xy (B « Bo) = ^frf^ (2-10) 

where we have invoked a small value of 7 to approximate the last numerator in Eq. ( 12.2ft by 
to 1. It is interesting to note that apart from the last factor which describes the entropy per 
coherence volume, this expression has the same parameter dependence as the one derived 
from Gaussian fluctuations in Ref. 0. 



III. MAGNETOHYDRODYNAMICS 

The remainder of this paper will revert to natural dimensionless units with h = k B = 
2e = v = 1, and absorb a factor of 1/c in the definition of B. 

Here we will focus on the nature of quantum critical transport in the hydrodynamic 
region 0, where fku <C ksl '■ The condition in Eq. ( 11.71) ensures Hlu c <C k B T, and a 
relativistic formulation is appropriate if also m < k B T is satisfied. We will use the method 
described by Landau and Lifshitz {46|, which was recently reviewed in the context of a 



string theory computation [47[. These previous analyses were carried out for B = and 
l/fimp = 0, and only considered the longitudinal electrical and thermal conductivities. Here 
we will show how the hydrodynamic analysis can be extended to include non-zero values 
of these parameters. Further, we will obtain results for the frequency dependence of the 
full set of transport co-efficients in Eq. ( II. 13ft . These results are consistent with the exact 
results obtained via AdS/CFT for a particular SCFT which are presented in Section IVl and 
Ref. E3— the latter results however extend over a wider regime of parameters. 

The fundamental ingredients of a hydrodynamic analysis are the conserved quantities and 
their equations of motion. Unlike in the theory of dynamics near classical, finite temperature 



critical points 49(, here we do not need to explicitly consider the order parameter dynamics 
for the effective equations of motion of the low frequency theory. The key difference is that 
ksT/h constitutes an intrinsic relaxation time for the order parameter fluctuations, and we 
are only interested in much lower frequency scales. In contrast, at a classical critical point, 
the relaxation time diverges. The frequency scales larger than k B T/h cannot be addressed 
by the methods below, and require a full quantum treatment of the dynamics of the CFT. 

The conserved quantities of interest are the electrical charge, the energy, and the mo- 
mentum. For the relativistic theories under consideration, these can be assembled into the 
electrical current 3-vector 2 J M = (p, J x , J y ), and the stress-energy tensor T^ v . We will use 



Notice that upon restoring v, the 3-current reads J M = (pv,0,0). 
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standard relativistic notation with the metric tensor g^ u = diag(— 1,1,1) and coordinates 
x^ = (vt,x,y). For the moment, we will ignore the momentum relaxation due to the weak 
impurity potential in Eq. (11.91) . and include its effects shortly below. With total momentum 
conserved, the equations of motion obeyed by the total electrical, momentum and energy 
currents are 



Q j(tot) M 



0. 



(3.1) 
(3.2) 



The first equation represents the conservation of charge and requires no further comment. 
The second equation describes the evolution of the stress energy tensor, and the term on 
the right hand side represent the effects of the external B field. Here F^ u is the applied 
magnetic field which takes the fixed value 3 



(3.4) 



and the right hand side of Eq. (13. 2p describes the Lorentz force exerted by this field, as 





(o 





o\ 










B 




v° 


-B 





discussed, e.g., in Ref. 
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In equilibrium, we have j( tot )^ = (p, 0,0) and then the term 
proportional to vanishes, as expected. 

To use the equations (13. ip and (13.21) . we need to relate j( tot )^ and T' tot ' F to parameters 
which define the local thermodynamic equilibrium, and a three- velocity which represents 
the ve 
usual 



ocity of the system in local equilibrium with respect to the laboratory frame. As 

-1, and = (1,0,0) in the 



461 ]. the three- velocity = dx^/dr satisfies u^u^ - 
equilibrium frame where there is no energy flow 4 . For the thermodynamic parameters we 



3 More precisely, if we also allow for an electrical field in the lab frame and restore v, the field tensor takes 
the form 

/ {c/v)E x (c/v)E y \ 
F» v = -{c/v)E x B . (3.3) 

\-(c/v)E v -B ) 

The equation of motion ()3.2|1 is a priori valid for velocities v' <C v, where the non-relativistic limit 
applies. However, it is valid as a relativistically covariant equation under Lorentz transformations 
characterized by the limit velocity v (not c), if it is understood that F^ v transforms as an antisymmetric 
tensor, F^ v = A^A^F"' 3 under a change of reference frame. This is exact up to negligible corrections of 
order 0([v/c] 2 ). 

4 Note that the absence of an energy flow in the rest frame defines the velocity vector u M in the " dynamic" 
frame used throughout this paper. Alternatively, one can formulate the hydrodynamics in the kinetic 
frame [45| where u M denotes the velocity associated with the electrical current via J M = pu^. However, 



the dynamic frame is the natural frame to work with, since we assume local equilibrium. u l is canonically 
conjugate to the momentum density T° l , and the energy current has a natural expression in terms of u M 
(cf. Eq. m» . 
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will use the charge density, p, the energy density e, the pressure P, and the magnetization 
density M; we define the pressure, P, as the negative of the grand potential per unit volume, 
and M as the derivative of the latter with respect to B. 

Using these parameters the stress energy tensor of a fluid is given by 



rp(tot) flL* 
rpflU 



_ ]1|C7 J?V _|_ rpEflU 

(e + P)u»u v + PgV + r^, 



where 



/o 





\ 
M 
\0 -M j 



(3.5) 



(3.6) 



is the polarization tensor 
of electric fields in the lab frame.) 
The electrical current is given by 



44| . (The electric polarizations M 0t = —M t0 vanish in the absence 



j(tot) M = jn + dvM ^ ; 



(3.7) 
(3.8) 



The 'dissipative current' accounts for the fact that the charge current and the energy 
current are not simply proportional to each other. This is because there is a heat flow even 
in the absence of matter convection, which is a consequence of particle-antiparticle creation 
and annihilation. 



We have introduced the transport currents 43j, J M and T^ u which represent observable 



quantities which can couple to probes externa 
j(tot)n ig the 

magnetization current 
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to the system. The remaining contribution to 



43j , which is induced due to spatial variations 



in the local magnetization density. The coupling of the magnetization to the magnetic field 
contributes an extra contribution —M^F" to the stress energy tensor, reducing its spatial 
diagonal to P mt = P—MB (see also Ap p endix IB]) . The origin and the physics of this term has 
also been discussed by Cooper et al 43j. Finally T Ellu represents the "energy magnetization 
current" . We will not need an explicit expression for this quantity here, apart from the fact 
that it obeys d v T E ^ v = 0. Expressions will be given later in the paper when we consider 
specific CFTs: for the super Yang Mills theory in Section |VJ and for the Wilson-Fisher 
fixed point in Appendix [Bl With these magnetization currents subtracted out, the residual 
transport currents continue to obey the equations of motion as in Eqs. (13.11) and (13.21) : 



d u T" u 



0. 



(3.9) 



In the expressions for the transport currents in Eqs. (I3.5P and (13.81) . we assume that e, P, 
p, and M are thermodynamic functions of the local chemical potential, /i, the temperature 
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T, and the magnetic field B. We will treat u^, //, and T as the "independent" degrees of 
freedom which respond to external perturbations, and assume that the remaining thermody- 
namic variables will follow according to the equation of state. In equilibrium, the non-zero 
components of and T^ v are 

J' = p ; T tt = e ; T xx = T vy = P — MB. (3.10) 



Eqs. (13. 5p and (13. 8p also contain t 
the current, as introduced in Ref. 
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re^dissipative components of the stress-energy tensor and 
these are are orthogonal to w M 

V "^«/ = 0, (3.11) 



and will be determined below by imposing the requirement that the total entropy increases 
under time evolution. 

We are now in a position to introduce the scattering due to a dilute concentration of 
impurities. We assume that their dominant effect is to introduce a relaxation in the local 
transport momentum density: impurity scattering conserves charge and energy, and we do 
not expect the magnetization currents to relax by impurity scattering. Thus, we modify 
Eqs. (E3D to 

V = °> 

dJT^ = F^J V + — (<f£ + vfuv) T^u T (3.12) 

'Hmp 

The new term in the second equation in Eq. (I3.12p represents the impurity scattering. The 
impurities, as described by the random potential in Eq. (11.91) are assumed to be at rest in 
the laboratory frame. The projection operators built out of the w M in the second term in 
Eq. (I3.2p ensure that in the laboratory frame only the total momentum, i.e., T l °, is relaxed. 
We will discuss a computation of the value of r imp later in Section HVl 

Following Landau and Lifshitz 46j, we now use the positivity of the entropy production 
to constrain the expression for the dissipative components and t^ u . First we notice from 
Eq. ([33]) that 

u u d^ = -(e + P)dX - u^e + u v d^\ (3.13) 
and from Eq. (13^12]) that 

u u d^ = F^u^ v . (3.14) 

Using the thermodynamic relations 

e + P = Ts + fip, de = Tds + fidp, (3.15) 

Eq. (13.1 ip and current conservation, Eqs. (13.131) and (I3.14p can be transformed into 

ii 1 t^ u 

a M (au") = - TfF^u^ - —d,u u , (3.16) 
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or 



0„ [su»-^ 



i {-d^ + F^u u ) + ^ 



T 



(3.17) 



It is natural to interpret the left hand side as the divergence of the entropy current. Ac- 
cordingly we can interpret the three vector 



sTu d 



(e + P)u^ - ^ = J E " - //J" 



(3.18) 



as the heat current. We have also introduced the energy current J ^ = (e + P)u^. 

Since the entropy can only increase, the right hand side of (13.1 7p must be positive. 
Generalizing the arguments of Landau and Lifshitz, we deduce the most general expressions 
for the dissipative currents which are linear in spatial gradients and the velocity, 



(-d uf i + F vX u x ) + /i^ 



-iiv 



(3.19) 
(3.20) 



Here rj and ( are the shear and bulk viscosities, and oq is a conductivity. Notice that 
there are only three independent transport co-efficients. We will neglect velocity gradients 
for the most part in this paper, and so the viscosities do not appear in our main results. 
Consequently, we have the remarkable feature that all transport response functions depend 
only upon a single dissipative transport co-efficient <jq. Notice that in the dissipative current 
the gradient of the chemical potential appears in combination with the electromagnetic forces 
F^Uv, which is natural since it is equivalent to an electric field. 



A. Linear response 



We will now follow the strategy of Kadanoff and Martin 50J: Use the equations of 
hydrodynamics to solve the initial value problem in linear response, and compare the results 
to those obtained from the Kubo formula in order to extract transport coefficients and their 
frequency dependence. 

First, we address the solution of the initial value problem in hydrodynamics. We begin 
by choosing our independent variables: from the structure of the above equations, it appears 
convenient to choose the four variables T, /i, and u x and u v . So we write 



/i(r, t) = fi + 5/i(r, t), 
T(r,t) = T + 5T(r,t), 

while F^ v is fixed at the value in Eq. (13.41) . We also write w M as 

/ 1 \ 



(3.21) 



v x (r,t) 
\v y {r,t) ) 



(3.22) 
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where v x , v y are of the same order as 5p and ST. 

The other variables, e, P, and p are constrained by local thermodynamic equilibrium to 
have the form 



p(r,t) = p + Sp = p + 

e(r, t) — e + 5e = e + 
P(r, t) = P + 5P = P + p5p + s5T. 



dp 


5p + 

T 


dp 


dp 


df 


de 


5p + 

T 


de 


dp 


df 



ST, 



ST 



(3.23) 



The various components of the stress-energy tensor and the current vector are perturbed 
accordingly. To linear order we have 



ST a = 


8e, 


ST U = 


rpti 


5T ij = 


5P5 


5J l = 


8p, 


5 J = 


J = 


V = 





-Vp + v x B + p 



VT 



(3.24) 



For small perturbations the conservation laws take the form 
d tP + V J = 0, 



d t e + V J 



0. 



dtP = -VP - r]V 2 v - ■ v) + J x B, 

with the energy and heat currents 
J E = (e + P)v, 

Q = J* 5 — pJ — (e + P)v — ppv — pv — sTv — pi>, 



(3.25) 
(3.26) 
(3.27) 



(3.28) 
(3.29) 



The crucial remnant of the relativistic theory in the linearized hydrodynamics is the fact 
that the energy and particle currents are in general not parallel, 



J* 



e + P 
P 

e + P 
P 



J — oq — V/i + v x B + p 



VT 



Tp 2 



-VP + Jx B) , 



(3.30) 



where we have used (13.151) to rewrite the dissipative current. The energy current consists 
of three parts: the first two are familiar from non-relativistic theory as the convection of 
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matter and heat flow due to a thermal gradient, with a thermal conductivity 46j k = 
(Tq{e + P) 2 /(Tp 2 ). The last term in ( [3.301) is proportional to the acceleration vector, and 
is a purely relativistic phenomenon |5l| . One can easily see that this term is responsible 
for the damping 7 of the cyclotron mode (cf., Eq. (11.191) ). by using the above relation to 
substitute for J in the momentum conservation law (13.271) . 

To complete the hydrodynamic analysis, we solve the equations (I3.25p - (I3.27I) for arbitrary 
initial values ST°, 5p° and = v® = 0, and obtain the response in the electrical current J 
and the heat current Q. The 'heat density' associated to the latter, q(r) = e(r) — pp{r) is 
canonically conjugate to the temperature at fixed chemical potential [50 1. 

After a Fourier transform in space and a Laplace transform in time, the linear response 
of any quantity A obeys 



A(k,u) 



G A;£ ^ p (k,uj)-G A;£ ^ p (k,0)5T°(k) t G A ; P (k,u)-G A ; P (W) 5 ofy 



iuj T 
G A . T oi(k,uj) — G A . T oi(k,0) 



IUJ 



)uJ 



(3.31) 



i=x,y 



where the coefficients are related to retarded equilibrium correlation functions, as can be 
shown from analyzing an adiabatic perturbation 50( of the form 

5T(r,t) 



5H(t) 



dr 



5n(r, t)n(r, t) 



T 



(e(r,t)-/in(r,t))-E^( r ' t ) T ° i ( r ' t ) 



.(3.32) 



Finally, using the conservation laws in the form 

iujp(k) = ikJ(k), 
iu{e{k) - fip(k)) = ik J^(fc) - fiJ(k) =ikQ(k), 

we obtain 



(3.33) 



A(k,u) = 



G A . (k,u)-G A . (k,O) ( VT°(k) 



1uJ 



+ G A]f (k,u;)-G A]f (k,0) ^ 



G A .T0i(k,L}) — G A -T0i(k,0) 



IUJ 



(3.34) 



i=x,y 



— * — * — # — * — » ^ 

For A = J and A = Q one recognizes the co-efficients of E = — V/i and (— VT /T) as 
(—1/iuj) times the Kubo formulae for the thermoelectric co-efficients a, 6c, k. The response 
to an initial velocity perturbation could be used to extract frequency dependent viscous 
response functions. 

After a Laplace transform in time, the Eqs. (13.251) and (I3.26P take the form 



dp 
dp 

de_ 
dp 



5p° + 



5p° + 



dp 
df 

de_ 
df 



5T C 



5T C 



UJ 



UJ 



dp 
dp 

de_ 
dp 



5p + 



5p + 



dp 
df 

de_ 
df 



ST) +iV (pv + Ba Q ev) - ia Q TV z 



ST + i(e + P)V ■ v 



2 (V 



(3.35) 
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for charge and energy conservation, and Eq. (13.271) 

i(e + P)v° = (u} + i/T imp )(e + P)v + i(pVp + sVT) 

- iBpev + iBa Q e(Vp - ^VT) + iB 2 a Q v 

+ ir]V 2 v + i(VCV-v), (3.36) 

for momentum conservation. 

In the case of weak enough momentum relaxation, the response functions will exhibit 
peaks associated with the normal modes of these linearized equations. Apart from the 
damped cyclotron mode discussed above, one finds two diffusive modes, as analyzed in 
Appendix [A] However, those will not be of importance below since we are restricting to 
long wavelengths in the sequel. 

In the limit k —>■ 0, the current and energy conservation impose 5p = i5p /u and 
5T = i5T°/u), expressing that the decay of initial perturbations can be neglected. Fur- 
ther, the contributions from viscosity can be neglected in this limit. Upon injection into 
the momentum conservation equations (I3.36p . we easily obtain the retarded Greens func- 
tions and via the mapping (I3.34j) and Kubo formulae the transport co-efficients defined in 

Eq. dnnj: 

(u + if r imp ) (lj + 27 + if r imp + ilj 2 / 7) 



&xx — ®Q- 



@xy 



OL, 



O^xy 



[(CJ + Z'7 -M/Vimp) 2 - 

p (7 2 + u 2 - 2iju + 27/r imp ) 



B [(lj + i>y + i/ r imp ) 2 - lj 2 ] 
s [lj c (ilj - l/r imp )(l - {^pp/{u 2 c Ts)){^ + l/r imp - iu))] 
B [(lj + r imp ) 2 - lj 2 ] 

s [ 7 2 + lj 2 + 7 (-^ + l/r imp )(l - pp/(Ts))} 
B [(lj + iy + i/r imp ) 2 - lj 2 ] 

e + P) 2 \ lj c1 



TBp J [(lj + i>y + i/ r imp ) 2 - lj 2 ] 



x 11 + (l/r lmp - M ^5|±2W) + (1/ _ ^ »V 



^xy 



7W*(e + P) 2 w imp ' lj 2 c (e + P) 2 f ' 

Ts 2 \ lj 2 



Bp ) [(lj + i>y + i/ r imp ) 2 - lj 2 ] 

i 2fia Q B(-f + l/T imp -iLj) f pa Q B \ 2 \ 
Ts lj c \ Ts ) J 

We also computed the thermoelectric co-efficients a xx , a xy by examining the heat current 
induced by an applied electric field, and precisely the same result as above was obtained. 
This confirms Onsager reciprocity which has to hold since the densities associated with 
the electric and heat currents are conjugate to 5p and 6T/T, respectively. The validity of 
Onsager reciprocity is a strong check of the consistency of our hydrodynamic description. 
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From expressions in Eq. (13.371) we can also derive the resistivities p = a 1 , the Nernst 
responses defined in Eq. ( ll.ISp and the thermal conductivities defined in Eq. (11 .151) . 

1 (uj + i/r imp )(uj + ij + i/r imp + iuj 2 /j) 

B {{ulhf + uj 2 - 2i{ulh)u + 2(c c 2 / 7 )/r imp ) 
p [(uj + iuj 2 /^ + i/r imp ) 2 -uj 2 ] 

s [(^/ 7 ) 2 + uj 2 + (uj 2 c h)(-iu + l/r imp )(l - p P /{Ts)) - (pp/(Ts)(-iu + l/r imp ) 2 ] ^ 
p [(uj + i-f + i/r imp ) 2 -uj 2 ] 

B (iuj-l/r imp ) 



Pxy 

u XX 

d 



xy ~ T [(uj + iulh + Vr imp ) 2 - uj 2 ] ' 
_ (e + P) (iu - l/r imp - 1^/7) 

T [(uj + iuj 2 /-f + i/r imp ) 2 -a; 2 ]' 

T [(uj + iulh + ,/r imp ) 2 - uj 2 ] ■ ^ 

These expressions contain the main results that were quoted in Section [H Although they 
appear rather complicated, most of the structure is tightly constrained and the predicted 
dependencies on to are robust consequences of hydrodynamics. 

Significant simplifications do appear if in addition to the small B assumption in Eq. ( 11.71) . 
we also assume that p is small. In particular, let us take B <C T 2 , p <ti T 2 and p ~ B ~ 
T ?j / 2 ^/uj. Note that in this limit we may simplify Eq. (13. 1 5[) to Ts = e + P. Then, the results 
in Eqs. (13.371) and (13.381) take the following more compact form (we have set l/r imp = 
because, as noted in Section fl] the dependence on impurity scattering is easily restored 
below by uj — > uj + z/r imp ): 

uj(uo + 27 + iuj 2 /^) p (7 2 + uj 2 — 2iqu) 

axx = aQ [(u + ijf-uj 2 ] ' **> = ~B [^ + z 7 ) 2 -^ 2 ]' 

1 u(u + iuj 2 c h + ij) B {{ulh) 2 + uj 2 - 2^/7)70;) 

,} -' a Q [{u> + iu>lh) 2 -ool] ' Pxy p [(uj + i uj 2 M 2 -uj 2 ] ' 

p iu s 7 2 + uj 2 — iju 



axx T [(uj + ry) 2 — uj 2 ] ' °^ B [(uj + z 7 ) 2 - uj 2 ] ' 

o s (uj 2 c h) 2 + uj 2 -i(uj 2 / 1 )uj _B_ iu 

XX P [(uj + iu 2 /j) 2 -u 2 ] ' xy T [(uj + iuj 2 /^) 2 — uj 2 ] 

_ iu — 7 _ uj c 

^xx ^77 '. 7n ITT 1 ^.7 



" M s [( w + ^ 2 /7) 2 -^ 2 ] ' Kxy "[(uj+iujihy-uji]- [6 - 6J) 

The above expressions are now easily observed to obey a remarkable 'self- duality' symmetry. 
Under the interchanges 

p <-> B and cr Q <-> 1/oq, (3.40) 
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the cyclotron frequency u c in Eq. f 1 1 . 1 8 j) remains invariant, while the damping frequencies 
in Eqs. (I1.19|ll.22l) interchange 

7 <->7u = ^c/7; ( 3 - 41 ) 
then note that the expressions for the transport co-efficients interchange as follows: 



@xxi @xyi C^xxi Qxyi ^xxi ^xy 



I 

Pxxi Pxyi $xy> ^xd ^xxi ^xy (3.42) 

These relations are consequences of the particle-vortex duality discussed in Ref. 0, and 
the mapping of the transport co-efficients in Eq. (I3.42p can be deduced from the mapping 
Ei <r-> €ijJj in Eq. (11.131) . These duality relations will be discussed further in the context of 



SCFTs solvable by AdS/CFT in Section [V] and in Ref. |3Cj: in this case the duality relations 
will be found to hold exactly for all p and B. 



IV. ESTIMATING THE MOMENTUM RELAXATION RATE 

As discussed in Section HJ we assume that momentum relaxation is caused by an external 
perturbation of the form 

Sim P = I dr I d d xV{x)0{x,r), (4.1) 



where V(x) is an external potential which is random function of spatial position x, but 
independent of r, with the averages in Eq. (11.101) . The operator O is the "thermal operator" 
of the CFT, i.e., the most relevant perturbation which drives the CFT massive (despite the 
name, it has nothing to do with temperature in the present context). 

We are interested in the modification of the equation of motion of the momentum density, 
T l ° by the impurity, because this is the only quantity whose conservation law is spoiled by 
the presence of impurities. In the absence of other perturbations from equilibrium we observe 
from Eq. (13. 1 2[) that the momentum density obeys 

°— = — - T %t + . . . , (4.2) 

uh ' imp 

where the ellipses indicates terms which have an explicit spatial gradient and so their spatial 
integral vanishes. We will describe here an estimate of r imp to order V^ mp . 

For definiteness, consider the Wilson-Fisher fixed point of a complex scalar ip in Eq. (11.31) . 
although the argument easily generalizes to other CFTs. We will also ignore the influence 
of B and /i as these are secondary perturbations. Then O = \tp\ 2 and 

T u = w*ditfj + c.c. , (4.3) 
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where w is canonical momentum conjugate to if). For the following, we need the commutator 

T' = i[2* 0] = a,(M 2 ). (4.4) 



We now compute l/ri mp using the memory function method 52|, |53j. From this approach, 
the estimate of the momentum relaxation rate is 

1 _ V-l P .. If d d k 



lim n - / T^M^-*, -^)T*(*, u)U , (4.5) 

where we are working in general d dimensions, and xt is the momentum density susceptibility 
i.e. 

nit ( _ _\/TTit 



XT = / d d xdT(T lt (x, r)T lt : (0,0)). (4.6) 
Using the scaling dimensions 

dim[T l °] = d + 1, 

dim[|^| 2 ] = d + 1 - 1/z/, (4.7) 

(where 1/z/ is the scaling dimension of the coupling conjugate to O = \ip\ 2 ) we obtain 

dim[x T ] = -d-1 + 2 dim[T i0 ] = d+l, 
dim [(T l {-k, -u)T\k, u))) ret ] = -d-1 + 2 dimpP] = d + 3 - 2/z/. (4.8) 

Thus ~ T d+1 , and 

1 ~K m p^ +1 - 2/ ^, (4.9) 



T-. lm P 

' imp 



which is the result quoted in Eq. (11.111) for d — 2. As noted there, v ~ 2/3, and so l/r imp 
is roughly temperature independent. 

An alternative, but less constructive argument proceeds along the lines of the discussion 
in Ref. |54|. From Eq. (14.11) we have dim[V] = 1/z/, and so from Eq. (jl.lOp we have dim[V| mp ] = 
— d + 2/u. This is indeed familiar from Harris' criterion which states that disorder is relevant 
if v < 2/d. Then, knowing dim[l/r imp ] = 1, the result (14. 9ft follows. 



V. DYONIC BLACK HOLES 

A. AdS4/CFT3 and the black hole solution 

From the point of view of studying quantum critical phenomena, the AdS/CFT corre- 



spondence [55| provides a wealth of new solvable strongly correlated conformal field theories 
(CFTs) in 2+1 dimensions. The key feature of these CFTs is that they admit a large N 
limit in which they can be described classically as a gravitational theory in 3 + 1 dimensions 
that asymptotes to Anti-de Sitter space (AdS^. The CFT is thought of as living on the 
'boundary' of the higher dimensional, or 'bulk', spacetime. 
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The correspondence furthermore allows us to consider departures from criticality due to 
finite temperature. This is precisely the type of systems we are studying in this paper. 
Finite tem per ature in field theory corresponds to allowing the bulk spacetime to contain a 
black hole 56|. The temperature of the field theory is just the Hawking temperature of the 
black hole. Finite temperature dissipation in field theory is dual to bulk matter fields falling 
into the black hole. 

We wish to consider a CFT with a global U(l) symmetry and a corresponding charge 
density p and a background magnetic field B. It was explained recently that this is dual 
to taking a dyonic black hole, carrying both electric and magnetic charge 29] . These black 
holes are solutions to Einstein-Maxwell theory in 3+1 dimensions. In this section, and also in 
a separate paper [30|], we will see how thermoelectric transport properties of the dyonic black 
hole precisely agree with our general analysis in the hydrodynamic limit. The black hole, via 
the AdS/CFT correspondence, provides a solvable example of the physics we are studying 
throughout this paper. Furthermore, the various particle-vortex dualities we have discussed 
above are seen to acquire a very transparent interpretation in the AdS/CFT correspondence. 

The canonical example of the AdS4/CFT 3 correspondence describes the infrared fixed 
point of maximally supersymmetric Yang-Mills theory with SU(N) gauge group. The dual 
gravitational theory in this case is M theory on AdS4 x S 7 [55(. However, in the large 
N limit and for the subset of field theory questions we are asking, this theory may be 
consistently truncated to Einstein-Maxwell theory with a negative cosmological constant in 
3+1 dimensions. The dimensional reduction is performed for instance in 14]. 

The action for Einstein-Maxwell theory with a negative cosmological constant may be 
taken to be 



V2N 3 / 2 
Gtt 



d 4 Xy 



-fj 



1 1 

4 4 



(5.1) 



which implies the equations of motion 



R 



2F„„F, 



fid*- V 



flU 1 



0. 



(5.2a) 
(5.2b) 



We have expressed the normalization of the action in terms of the field theory quantity N. 
A black hole in AdS& with planar horizon has metric 

C( 2 r „, s ,1 , o , on 1 dz 2 

J 2 



ds 2 = ?L [-f(z)dt 2 + dx 2 + dy 2 ] 



z 2 f(z) ' 

The dyonic black hole carries both electric and magnetic charge 

F = ha 2 dx A dy + qadz A dt , 



(5.3) 



(5.4) 



where h, q and a are constants that will be related to field theory quantities shortly. The 
function appearing in the metric is then 



f(z) = l + (h 2 + q 2 )z 4 -(l + h 2 + q 2 )z 3 



(5.5) 
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In these coordinates, the conformal boundary of the spacetime is at z — > whereas the black 
hole event horizon is at z — 1. 

In the following subsection we summarize the thermodynamic properties of this black hole 
spacetime that were derived in which is also the thermodynamics of the CFT. We will 
express thermodynamic quantities in terms of the dual field theory background magnetic 
field and chemical potential. These are related to the constants q and h that appeared in 
the black hole solution as follows 29 

B = ha 2 , n = —qa . (5.6) 



B. The grand canonical ensemble 



We give the thermodynamics of the CFT in terms of the temperature T, the chemical 
potential p and the background magnetic field B. Many variables are most simply expressed 
in terms of an auxiliary quantity a(T,p,B) which is determined from 



AttT 



a 



,,2 r>2 

3 _ H £_ 

a 2 a 4 



The thermodynamic potential is 

67r 4 \ a" a 
Here V = fdx dy is the spatial volume. The energy density is 



a 2 B 2 
-1-^ + 3— 



V2N 3 / 2 a 3 

6v y 



, P 2 B 2 

a 2 a 4 



The entropy density is 
The charge density is 
The magnetization is 
The pressure is 



6~ 



-a 



M 



y/2N 3 / 2 2f i 

p = a — . 

67r a 



1 on V2N 3 / 2 b 

-a— 

or 



VdB 



671 



P = MB + - . 

2 

These formulae satisfy the thermodynamic relation 

e + p = Ts + y.p • 



(5.7) 



(5.8) 



(5.9) 

(5.10) 

(5.11) 

(5.12) 
(5.13) 

(5.14) 
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Note that the above implies the relation 



Finally, it is useful to define 
Which gives the relation 



£ = -(sT + np-MB) . (5.15) 
3 



X=-toT-T. (5.16) 



P MT 2 = -fiB X 2 . (5.17) 



C. Magnetization densities and currents 

In the following subsection, we will use Kubo formulae to the obtain transport coefficients 
of the SCFT from retarded Greens functions. When applying the Kubo formula to systems 
with a background magnetic field, it is necessary to subtract effects due to magnetization 



currents, as explained in Cooper et al. [43J]. The magnetization currents are 

JT* = *jf>jM , (5.18) 

and 

T^ g - = e ij d j M E . (5.19) 

Here and below the indices i, j run over the two spatial coordinates x and y. The equilibrium 
magnetization density M and energy magnetization density M E for the dyonic black hole 
are obtained by differentiating the free energy with respect to a constant magnetic field for 
either the charge or momentum currents, 

M = ~, (5.20) 

u± xy 

M =-Jf^- (5 - 21) 

xy 

Here we define SF E y = d x 5g^ y — d y 5g^ x and 5F xy = d x 6Ay — d y 6A®, where 5g® a is a background 
gauge field sourcing T ta , and 5A° a sources J a . Further comments on these magnetization 
densities, and their computation for the scalar field theory in Eq. (11. 3ft appear in Appendix iBl 

In the AdS/CFT correspondence, the free energy Q is just the action of the dual gravita- 
tional background. To compute the derivatives in (15.201) and (15.211) we must consider on shell 
fluctuations of the bulk metric and gauge field that tend towards 8g® a and 6 A®, respectively, 
near the conformal boundary z —* 0. We then differentiate the action with respect to the 
boundary values of these fields. 

More concretely, the boundary condition may be taken to be 

5 A y -> xB as z^O, (5.22) 
5g ty z 2 /a = 6Gy — > xB E as z , (5.23) 
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with B and B E constants and all other fields having normalizable behavior near the bound- 
ary. It turns out that the linearized fluctuation equations about the dyonic black hole 
background with these boundary conditions may be consistently truncated to the fields 8A y , 
8 At and 8gt y . The solution is 

8 A y = x(B - qB E z) , (5.24) 

SG y = xf(z)B E , (5.25) 

_ , hB E , o . hB . , 

8A t = —(z 2 - 1) (z - 1) . 5.26 

la qa 

Note that A t vanishes at the horizon, z — 1, as required. 

Because we are fluctuating about a solution, the linear variation of the bulk action will 
vanish. However, there will be a boundary term that arises due to integration by parts when 
evaluating the action on the solution. There is also a boundary term that must be included 
to renormalize the gravitational action 



bdy. 



y/2iV 3 / 2 
Qtt 



(5.27) 



where 7 is the boundary metric and 9 = r y ,iV 6 ia/ is the trace of the extrinsic curvature 
9^ u = —^(V^n u + V^n^), with n an outward directed unit normal vector to the boundary. 
For fluctuations about a solution, we have 



81 = — — — / d d xV-7 



-F^SA" + \ (9^ - 9r v ~ 27^) 8g»v 



(5.28) 



Evaluated on the background (I5.3P and (15.41) . and considering only 8A y , 8A t and 8g ty , this 
expression becomes 

v/2iV3/ 2 , f. 



81 = qa z / d A x8A t . (5.29) 

6tt J 

Note that only the first term in the variation of the action (15.281) contributes. Using the 
solution (15.261) for 8A t , we obtain the magnetization 

8B 6n 6ir a 

in complete agreement with our previous expression (I5.12p . For the energy magnetization 
we obtain 

. r w 8S qa , , uM 

M = —a- — - = — — M = — — . 5.31 
8B E 2 2 v ' 

D. Transport coefficients in the d.c. limit 

We will obtain the transport coefficients using Kubo formulae for the retarded Greens 
functions. The Greens functions are obtained by considering fluctuations about the dyonic 
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black hole background. In 29] these functions were obtained at k = and to leading order 
as u — > with B and p held fixed. Unlike in our general magnetohydrodynamic (MHD) 
analysis in Section IHIl no assumptions are made here requiring B to be small. The current- 
current correlator is 

= -;§*«/«&, (5-32) 

the current-momentum correlator is 

3e 
2B 

and the momentum-momentum correlator is 



G% t] (uj) = -—iueij , (5.33) 



xT 3 s 2 9pe 2 T 2 

G T u T tj M = ~ x2B2 + p2T2 ^, - AB{x 2 B 2 + p 2 T ^ ■ ( 5 - 34 ) 

To relate these results to our general MHD study, we need an expression for the heat 
current Q l . This is defined in Eq. ( I3.18P ; using the expression for the stress energy tensor 
in Eq. (13.241) . we see that in linear response (small velocities with respect to the lab frame) 
we can work with 

Q i = T u - fiJ* , (5.35) 

leading to the correlators 

<W*M = (^- + m)uj€ 1Ji (5.36) 



and 



rR f v XT 3 S 2 . -ps 2 T* + B 2 p 2 p X 2 + pT 2 M 2 B 2 , 

W<") = - x 2 B 2 + p 2 T 2 ^ + B( X 2 B 2 + p 2 T 2 ) • (5 ' 37) 

In obtaining these expressions, we used (15.151) and (15.171) . 
The electrical conductivity is given by the Kubo formula 

ImG^ Ti (uj) p , 

a = - lim £f±l = !L e . (s.gg) 

The other thermoelectric tensors are also given by a Kubo formula. However, we should 

use the transport currents which are obtained from the supergravity currents by subtracting 
the magnetization currents [43] . The correct Kubo formula is 

1 lmG% Q3 (uj) M 

a y = -- hm + . (5.39) 

Thus we obtain 

a H = gCij ■ (5-40) 
Similarly, the heat conductivity is given by the Kubo formula 



1 ImGg^w) | 2(M E -pM) 
T ui 



^ = -±b lI2E + 2^ja % , (5.41) 
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or 

xT 2 s 2 ps 2 T 4 
^ = X 2 B 2 + f?T* Sij + TB( X 2 B 2 + p 2 T 2 f J ' (5 ' 42) 
We can now compare these results with those of our general MHD computations, by 
taking the uj — > limit of the MHD transport coefficients (I3.37p . We see immediately that 
the expressions for aij and a^- agree exactly. In order to match 7^- we need to recall that the 
MHD results are only valid for small magnetic fields B ^ T 2 , see Eq. ( II. 7p . Furthermore, 



we need to know the conductivity oq for the dyonic black hole. It is shown in Ref. [30 that 
for the black hole 

s \ X 



<">\**=\j+p)%- (5 ' 43) 

Using this formula and taking the small magnetic field limit, we obtain 

a Q {e + P) 2 s 2 T 
K *il_B<T 2 ~~ v ~B~p €ij ' (5.44) 

This expression now agrees exactly with the corresponding limit of the MHD result (13.371) . 
Thus we see that the dyonic black hole fits into the general class of finite temperature 
deformations of quantum critical points that we have studied via hydrodynamics. The black 
hole gives specific values for <jq and the other thermodynamic quantities and furthermore 
allows the results to be extended to arbitrary magnetic field. 



E. Bulk electromagnetic duality and CFT particle/ vortex duality 

A consequence of the bulk description is that it gives a very transparent rationale behind 
the dualities in the transport coefficients that we commented upon earlier. The study of 



electromagnetic duality in the dyonic black hole is pursued in depth in Ref. [30j, which 
furthermore obtains expressions for the black hole transport coefficients away from the d.c. 
limit. Here we shall summarize some results from that paper and show how they precisely 
match our expectations from MHD. 

The bulk Maxwell theory enjoys an electromagnetic duality. This interchanges the bulk 
electric and magnetic fields E — > B and B — ► — E. Acting on the dyonic black hole solutions 
(15. 3p . this corresponds to q — > h and h — > —q. Using (I5.6P and (15.111) and the fact that the 
bulk coupling is inverted under electromagnetic duality, this implies that 

Tp TB X T 
B^ — , , =->-. 5.45 

X X T x 

Thus we see that the bulk electromagnetic duality maps the CFT into the same CFT with 
the values of the background magnetic field and charge density interchanged. This is the 
origin of the particle-vortex duality that we noted in our MHD computations. Indeed it is 
immediately seen that under the transformation 1 15. 45ft our expressions for a, a and R in Eqs. 
(I5.38p . 1I5.40P and 115.421) transform according to 113.421) . There are some overall factors ofx/T 
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different to (13.421) . due to the fact that the transformation (I5.45P is normalized differently 
to f)3.40p . The normalization in (15.451) is natural from the string perspective. 

It remains to see how the thermoelectric transport coefficients of the black hole transform 



under this map. This is shown in detail in Ref. [3CJ. The central point is the following. 
The bulk Maxwell potential A determines the bulk electric and magnetic fields through 
E ~ d t A and B ~ d z A. As we tend towards the boundary z — > 0, the electric piece is 
non-normalizable, and results in a boundary electric field E. The magnetic piece however 
is normalizable and therefore results in a boundary current J. The bulk electromagnetic 
duality is thus seen to exchange the electric field in the CFT with the current. More precisely, 
one finds 

Ei <-> eijJj . (5.46) 

As we commented below Eq. (I3.42p . this map together with the definition of the transport 
coefficients in Eq. (11.131) is enough to obtain all the duality transformations (13.421) . 

The results for the transport coefficients presented in the previous subsection were only 



obtained in the d.c. limit uo — > with B and p fixed, following Ref. |29|. However, using 
the AdS/CFT dictionary, it is possible to study thermoelectric transport at all frequencies. 
This is done in Ref. |30|. In particular, taking the limit uj — ► with p ~ B ~ T 3 ^ 2 y/uj, one 
obtains precisely the same expressions as those that followed from the MHD analysis (13.391) . 
thus providing a consistency check for our picture. One can go further with the dyonic black 
hole and study transport and all p, B and u numerically. For all values, the particle-vortex 
duality holds automatically because of the bulk electromagnetic duality. This is the power 
of the AdS/CFT correspondence: all transport phenomena of the strongly correlated CFT 
at large N are reduced to solving the equations for classical perturbations of the dual black 
hole in Einstein-Maxwell theory. 



VI. CONCLUSIONS 

This paper has presented a general theory for hydrodynamic thermal and electric trans- 
port in in the vicinity of a quantum critical points described by "relativistic" quantum field 
theories. We have also shown that the results constitute a valuable starting point to under- 
stand experimental observations in a regime where no previous description was available. 

It is perhaps useful to describe the results here from the vantage point of the Galilean- 
invariant hydrodynamic approaches which are traditionally used in condensed matter physics 



50, [57|. In such theories, charge (or number) currents are proportional to the momentum 
current, and consequently the conductivity is infinite in the absence of impurities (in the 
presence of a magnetic field, this implies Kohn's theorem js^]). The natural transport co- 
efficient is the thermal conductivity, and this determines various diffusivities and damping 
constants. 

In contrast, in the present paper, we have used a very different starting point. We 
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considered a theory with both particle and anti-particle (hole) excitations, in which there is 
no proportionality between momentum and charge currents. For the case with particle-hole 
symmetry (p = 0) and B = 0, the momentum and charge currents are decoupled from each 
other, and it is possible to have a charge current with no momentum current: the electrical 
current can decay to zero from such a state, and this decay is associated with the universal 
electrical conductivity <tq 17|. There is no analog to <jq in the Galilean invariant case. Also, 
because of the symmetry of the stress-energy tensor, we could identify the energy current 
with the momentum density; the conservation of momentum then implied that the thermal 
conductivity was infinite 58(. Upon relaxing the requirement of particle-hole symmetry {i.e. 
allowing p ^ 0), we found the appearance of some characteristics of the Galilean invariant 
systems. In Eq. ( 11.21) we found that the excess particles (or holes) contributed a Drude-like 
conductivity above the quantum-critical <jq. The thermal conductivity became finite, but 
with a value related to oq by a Wiedemann-Franz-like relation. Finally, we also turned on 
a B ^ 0, and showed how all of the longitudinal and transverse transport co-efficients could 
be related to <tq in relationships that were summarized in Section [H 

While our analysis was specialized to relativistic quantum critical points, we expect that 
many aspects will generalize to other strongly interacting quantum critical points. Only a 
discrete particle-hole symmetry is required to decouple the charge and momentum currents, 
and this should be sufficient to obtain a finite oq. 

We also discussed applications of this general hydrodynamic structure to measurements 
of the Nernst co-efficient in the cuprates and NbSi films. For reasonable values of the 
parameters, we were able to reproduce several key aspects of the B and T dependence 
of the observations. Our results also make a significant prediction, characteristic of such 
"relativistic" theories of the superfluid-insulator transition: the presence of a hydrodynamic 
cyclotron mode. For the simplest case of a superfluid-insulator transition of Cooper pairs 
at integer filling as described by Eq. (11.31) . this cyclotron mode can be considered due 
to the motion of Cooper pairs (or their dual vortices). However, for the more complicated 
examples at fractional filling noted in Section [TAJ such a simple interpretation is not possible, 
and the cyclotron mode is due to motion of all charge carriers, including those carrying 
fractional charges. From our fits to the data in the cuprates in Section [III we found that 
in presently studied samples this cyclotron mode was strongly overdamped by impurity 
scattering. However, this raises the possibility that the cyclotron resonance which might be 
observable in ultrapure samples. We estimated that the hydrodynamic cyclotron frequency 
in the cuprates was smaller than the cyclotron frequency of free electrons by a factor of 
order 10~ 3 . Observation of this resonance would constitute a strong test of the theoretical 
ideas presented here, and we hope such experiments will be undertaken. 

Another class of results in this paper described the remarkable holographic connection 
between this intricate hydrodynamic behavior in 2+1 dimensions and the quantum theory 
of dyonic black holes in 3+1 dimensions. Using the AdS/CFT connection, we presented in 
Section [V] exact results for the hydrodynamic response functions of the vicinity of a 2+1 
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dimensional supersymmetric conformal field theory. In the appropriate limiting regime, 



these results, and those in Ref. |30|, were in complete agreement with those obtained from 
the magnetohydrodynamic analysis in Section IIHI This agreement lends strong support to 
the validity of our MHD analysis. Indeed, the analysis of the dual gravity theory helped 
guide our determination of the MHD equations. 
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APPENDIX A: NORMAL MODES OF THE HYDRODYNAMIC EQUATIONS 



It is interesting to analyze the normal modes of the linearized magnetohydrodynamic 
equations in the absence of impurity scattering. Assuming a space and time dependence 
e -iut+ik-x Q f ^ fiji an( j ._ m anc j y ^ ._ (fclty . (eu), we find four independent modes 

satisfying the equations 



de_ 
dp 



. de 



5T - k(e + P)v\\ 







dp 
dp 



dp 



6T^j - k (pun + Ba QVl _) + ik 2 a Q (dp - ^5T) = , 

oo(e + P)v\\ - k(p5p + s5T) - ipBv ± + iB 2 a Q v\\ - ik 2 (r] + ()v\\ = , 
u(e + P)v± + kBa Q (dp - ^5t) + ipBv\\ + iB 2 a Q v± - ik 2 r] Vl _ = . (Al) 

In the long wavelength limit fe<l, one finds two modes corresponding to damped cyclotron 
oscillations of the plasma 

u± = ±u c — ry . (A2) 

These modes have a velocity field with v\\ = ±iv±, while 8p and 5T are smaller than v\\,v± 
by a factor of order 0{k). 

Further there is a diffusive mode with frequency proportional to the conductivity <jq and 
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a quadratic dispersion relation 



k 2 a Q (e + Pf 



( de 


dp 
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dp 




[ail 
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(A3) 



T 



<9e 

This mode has no fluctuations in energy density, 5e = — — 

dp 



(P 



2 ' B 2 o 2 



. de 
6 » + df 



5T = 0. The velocity 



field v is of order 0{k) relative to dp, 5T. 

Finally there is a subdiffusive, transverse shear mode with strongly suppressed fluctu- 
ations in temperature and longitudinal velocity component 5T = 0(k 2 ), v» = 0(k 3 ). It 
exhibits an unusual dispersion relation 



^subdiff 



lk 4 7] 



dp 
dp 



(A4) 



and we have the relation 



ikdp w Bv±. (A5) 

The origin of the k 4 dispersion ( 1A4I) can be seen as follows: A strongly suppressed u(k) 
implies that momentum density is nearly conserved. Hence the total force density vanishes 
to lowest order, i.e., the Lorentz force is balanced by a longitudinal pressure gradient and a 
transverse friction force, 



(Jx B)\\ w VP + C%C,w), 
T]\7 2 v ± m (Jx B)± = -BJn. 



(A6) 
(A7) 



The first equation yields relation (1A5j) . The second can be injected into the equation for 
charge conservation 



d t p 



-V • J w V 



5 



V 2 K) 



^ w2 / V/i 



5 

from which the dispersion follows upon using d t p = —iu — 

dp 



-y 2 



B 



(A8) 



5p. 



APPENDIX B: MAGNETIZATION AND ENERGY MAGNETIZATION 

In our computation of the transport co-efficients using the Kubo formula in Section |V] 
we had to face the issue of the subtraction of magnetization currents, as discussed earlier in 
Refs. 
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These subtractions were computed in Section [V] using the AdS/CFT mapping. 
This appendix describes the nature of these magnetization subtractions in the context of the 
scalar field theory in Eq. (jl.3p . Actually, most of the basic issues are already clarified in free 
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field theory, and so we will limit our discussion here to this simple case. The generalization 
of the free field results to the interacting Wilson-Fisher fixed points can be straightforwardly 
carried out along the lines of Refs. 



31 



32, and so we will not discuss it here. 



So we consider here the free field version of Eq. 11.31 with Lagrangian 

C = [(d» + iA»)r\ [(fy - iAM + m 2 \^\ 2 . 



The stress-energy tensor is 59] 



T, 



[(<9 M + iApW] [{d v - iAM + \{d u + iA u )r\ [{d„ - iA, 



while the U(l) current is 



J M = (dp - iA^ip - i^dfj, + iA^ip* 



The equation of motion is 



{d^-iA^-iA^^m 2 ^. 



(Bl) 



(B2) 



(B3) 



(B4) 



It is now a straightforward, but tedious, exercise to verify that the above expressions do 
indeed imply the MHD equation of motion in Eqs. (13. 1113.21) . 

For the thermodynamics, we need the particle and hole eigenenergies. These are organized 
in Landau levels, with energy 



e 2 i = 2B{£+l/2) + m 



(B5) 



(£ = 0,1, . . . oo) and degeneracy per unit area of B/(2ir). From this, we can easily obtain 
expressions for the grand potential (for all thermodynamic quantities we subtract an infinite 
T = value, and u n is a Matsubara frequency which is an integer multiple of 2n): 
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BT 



2tt 
BT 



— Y [ln(l-e- 
2vr 1 v 



(«/-f*)/T 



) + In (1 



)]■ 



(B6) 



We also obtain the entropy, s, the density, p, and the magnetization density, M, by 

i an 



S= VdT ' 9 



idn i an 

' ~~VdB 



(B7) 



Following Cooper et al. |43j, it is useful to define an internal pressure Pj nt which equals 
P int = P-MB 

B 2 ^e + i/: 



+ 



1 



B8) 



39 



We define the energy density e by (T a ). Evaluating this from (1B2|) in Euclidean Matsubara 
space, we obtain 



(Ttt) 



2tt 

B 
2^ 



E 



-(u n - ifi) 2 + e\ 
(uj n -ifj,) 2 + ej 



+ 



1 



(B9) 



It is now easily verified that the relations in Eq. (I3.15P are obeyed. 
In a similar manner, we can compute (T xx ) and find 

ifi) 2 — m\ 



(T xx ) 
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Pi 



int; 



(BIO) 



and so Eqs. 113. lOR are also obeyed. 

Let us now write down the explicit result for M from Eq. (IB7I) : 



M 



(ei-A*)/T 



5 i+ 1/2 



2^ 



) + In (1 - e 
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e (e t -n)/T _ X e (e e +li)/T _ X 



(Bll) 



A direct evaluation of (r x J) in an infinite system yields only the second term but not the 
first. We now argue that the first term is the contribution of edge states. Notice that this 
first term can be rewritten in the form 



M 



2tt J 



dEg e (E) 



B ^ £+ 1/2 



e (E-n)/T _ x 
1 
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t 1 

where we can interpret g e (E) as the magnetization of edge states: 
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9e(E) 



for E < e 
for €£-\ < E < ee 



(B13) 



Note that g e {E) is a piecewise constant function, and there is one additional edge state as 
each Landau level is crossed, as expected. We can write the expression for the magnetization 
as a sum of a bulk and edge contributions as 

1 



1 f°° 
M= 2^J dEg{E) 



+ 



e (E+n)/T 



(B14) 
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where 



9 (S) = «.(£)- f) S£±^(£--„, 



£=0 



(B15) 



With the above form for M, we can now immediately use the results of Ref. |42j (compare 
their Eqns. (31) and (34)) to obtain the value of M E : 



— I dEEg(E) 
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e (E-n)/T _ I e (E+ti)/T _ i 



(B16) 



The subtraction for n xy is, by Cooper et al. 43J Eq. (69), 2M® where 



M Q = m e - /iM 
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ee+fl EdE 
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.(B17) 



An alternative evaluation of M®, which is generalizable to interacting theories, follows 
from the representation of the heat current discussed in Ref. |6C|. As noted in Section IV CI we 
need the response to a "magnetic field" which associated with the energy (or heat) current, 
just as the ordinary magnetization is a response to a magnetic field associated with the 
charge current. So we introduce a vector potential Aq which couples to heat current: this 
is done by the replacement [60(| A — > A + iAqUu. Consequently, the free energy density in 
the presence of this additional "magnetic field" Bq is obtained from Eq. flB6|) simply by the 
replacement B — > B + iBnu n , which yields 



n t 



y = — J2( B + iB Q^n) In [K - tfi) 2 + 2{B + iB Q u n )(£ + 1/2) + m 2 ] 



V 2tt 

Taking the derivative with respect to we obtain 
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which agrees with Eq. (1B17|) . 
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